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SECTION  I 


INTRODUCTION 

The  earliest  model  proposed  for  the  determination  of  cumulative  fatigue 
damage  gave  a  result  which  is  now  called  Miner's  rule  (see  [1])  despite  the 
presentation  of  this  formula  for  life  prediction  20  years  ea  lier,  in  [2], 
by  A.  Palmgren.  It  is  a  deterministic  model  based  on  a  linear  damage  theory 
which,  strictly  interpreted,  has  been  known  to  be  false  for  a  long  time.  Its 
acceptance  as  an  engineering  tcol  in  design  has  been  hard  to  overcome  because 
of  its  ease  of  application  and  the  many  situations  in  which  the  predictions 
are  reasonably  good.  In  fact,  it  was  the  best  qualified  of  all  those  that 
were  compared  at  one  time;  see  references  in  [3]. 

The  generalization  from  that  simple  deterministic  model  has  proceeded 
mainly  on  two  separate  fronts.  In  one,  it  has  been  replaced  by  a  probabilistic 
framework  in  which  fatigue  life  is  regarded  as  being  stochastic  in  nature.  In 
these  models  the  statistical  problems  of  determining  the  family  of  distributions 
and  the  regression  of  mean  stress  on  life,  along  with  the  estimation  of  para¬ 
meters,  have  come  to  the  fore.  Only  recently,  from  such  a  model  for  fatigue 
crack  growth,  has  Miner's  rule  been  shown  to  be  the  expectation  of  fatigue  life 
(which  is  stochastic)  under  repea' ed  load  spectra. 

Under  more  realistic  conditions  than  those  assumed  by  Miner,  namely 
that  cumulative  damage  is  dependent  on  load  order,  a  modified  Miner's  formula 
was  obtained  which  was  theoretically  capable  of  taking  into  account  stress 
history.  This  result  yields  a  formula  for  the  expected  number  of  spectra 
which  can  be  repeated  until  failure.  It  contains  as  a  special  case  the  orig¬ 
inal  Miner'r  rule  with  interactions  neglected.  A  series  of  publications  which 
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give  this  development  are  [4],  [  5  1.  [  3  ]•  The  article  [  3  ]  is  of  an  exposi¬ 
tory  nature,  providing  a  broader  scope  with  some  historical  perspective. 

On  the  other  front  there  have  been  a  great  number  of  investigators 
studying  the  physical  mechanisms  of  crack  growth  using  various  mechanisms  in¬ 
volving  fracture  and  dislocation  theory  of  metals.  Some  of  these  which  explore 
in  detail  the  consequences  of  load  order  interactions  are  [6],  [ T  J »  [8],  [9]. 

One  of  the  latest  is  Impellizzeri  [7]-  That  paper  gives,  albeit  in 
a  deterministic  form,  a  cumulative  damage  rule  which  can  be  used  to  predict 
the  number  of  spectra  which  can  be  repeated  until  failure,  taking  into  account 
the  load  interaction  factors.  It  is  a  modification  of  the  Miner's  rule  sum 
and  agrees,  except  for  notational  differences,  with  the  general  results  in 
[4],  [  5  ]  for  the  expected  value  of  a  stochastic  crack  growth  model. 

Most  of  the  present  investigators  assessing  cumulative  fatigue  damage 
emphasize  only  one  or  the  other  of  1)  identification  of  the  important  de¬ 
terministic  factors  in  crack  initiation  or  2)  statistical  problems  associated 
with  the  pooling  of  the  variation  within  material,  manufacturing,  loading,  and 
environment  into  the  several  unknown  parameters  of  a  statistical  law.  Neither 
of  these  by  itself,  we  feel,  will  ultimately  be  adequate  for  a  complete  theory. 
What  is  now  needed  is  a  more  comprehensive  interdisciplinary  approach.  The 
statistical  theorist  must  utilize  the  knowledge  gained  by  the  investigators  who 
have  concerned  themselves  with  quantifying  the  effects  of  such  factors  as  stress 
intensity,  residual  stress  and  the  interaction  of  varying  peak  stress  levels 
in  a  loading  sequence;  see  [7]  and  the  reference  there.  Likewise  the  fracture 
mechanist  must  also  utilize  the  knowledge  which  has  been  gained  from  the  mathe¬ 
matical  and  statistical  methodjlogy  of  fitting  parameters  to  life  distributions. 
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That  this  problem  is  well  recognized  can  be  seen  from  the  comments  of 
J.  Schijve  which  were  made  during  the  second  F.  J.  Plantema  Memorial  Lecture, 
given  at  the  11th  conference  of  the  International  Committee  on  Aeronautical 
Fatigue  held  in  Stockholm  in  May  1969  (see  [10]): 

"Damage  theories  so  far  have  been  unable  to  include  interactions 
between  fatigue  loads  of  different  magnitudes  (sequence  effect)  whereas 
laboratory  tests  have  clearly  indicated  that  such  interactions  do  occur 
and  can  be  qualitatively  understood." 
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SECTION  II 


THE  MECHANISMS  OF  FATIGUE  CRACKING 

Before  an  accurate  theoretical  model  of  the  fatigue  mechanism  is 
constructed  it  is  well  to  summarize  some  of  the  empirical  results  which 
have  been  obtained  by  earlier  investigators.  Thus  we  hope  to  make  certain  that 
the  abstract  model  agrees  with  particular  facts.  In  order  to  present  a  complete 
picture  of  the  state  of  knowledge  we  paraphrase  from  the  appropriate  sections 
in  the  more  inclusive  report  of  McMillan  and  Pelloux  in  [11],  rather  than  [8]: 

The  influence  of  maximum  stress,  stress  range,  and  sequence  of  load 
application  on  the  rate  and  mechanism  of  fatigue  crack  propagation  in  2024-T3 
aluminum  alloy  was  carefully  studied  by  means  of  electron  fractography. 

Variable  amplitude  loading  programs  were  designed  to  provide  tests  under  the 
following  conditions: 

1^  Constant  maximum  stress  with  three  different  levels 
of  stress  range, 

2)  Constant  stress  range  with  three  and  four  levels  of 
maximum  stress, 

3N  Pseudorandom  load  application  achieved  by  random  dis¬ 
tribution  of  the  load  spectra  defined  in  items  l')  and 
2). 

U'l  Uniform  maximum  stress  with  peak  overloads  and  under¬ 
loads. 

The  macroscopic  growth  rates  were  determined  on  center-notched  crack 
growth  panels  and  the  fracture  surfaces  were  examined  by  electron  fractography. 
The  analyses  of  the  influence  of  program  loads  on  the  rate  and  mechanism  of 
fatigue  crack  growth  were  accomplished  by: 
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1)  Comparing  plots  of  crack  length  versus  measured 
rates  of  crack  propagation  for  the  different  pro¬ 
grams.  The  measured  rates  were  also  compared  with 
rates  calculated  by  a  computer  program. 

2^  Relating  the  count  and  spacing  of  the  fatigue  stri- 
ations  observed  on  the  fracture  surfaces  to  the 
applied  load  program,  by  means  of  electron  frac- 
tography . 

These  carefully  controlled  fatigue  crack  propagation  tests, made 
under  program  and  random  loads,  were  followed  by  a  study  of  striation  cowits 
and  spacing  measurements  using  electron  microscopy,  all  of  which  help  to  pro¬ 
vide  a  better  understanding  of  the  mechanisms  of  fatigue  cracking.  The  main 
findings  with  respect  to  aluminum  alloy  2024-T3  are  as  follows: 

1)  The  advance  of  the  fatigue  crack  front  takes  place 
only  during  the  stress  rise  portion  of  a  cycle. 

2)  The  profile  or  sides  of  a  striation  are  related  to 
the  loading  and  unloading  sequence  of  a  load  cycle. 

3)  Microscopic  crack  growth  rates  were  measured  as  a 
function  of  load  amplitudes,  AS  ,  and  maximum  loads, 

S  max,  and  an  empirical  crack  growth  rate  equation 
was  obtained: 

crack  growth  rate  constant  (AS)n(S  max)m  . 

This  equation  relates  the  relative  growth  rates  at 
a  given  crack  length  under  program  loading. 
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4)  In  random  load  programs  with  constant  maximum 
stress,  the  growth  rate  is  the  same  as  the  pro¬ 
grammed  spectra  rates.  No  sequencing  effect  is 
apparent  in  changes  of  stress  variation  alone. 

5'  In  random  load  programs  with  variable  maximum 
stress,  the  crack  growth  rates  were  higher  than 
for  the  equivalent  programmed  load  programs.  The 
importance  of  the  crack  tip  radius  in  controlling 
crack  rates  for  random  loads  was  pointed  out. 

6)  The  measured  crack  growth  rate  of  program  and  ran¬ 
dom  load  crack  propagation  tests  were  compared 
with  rates  calculated  by  a  computer  program.  The 
computed  rates  are  not  conservative  but  they  do 
not  differ  from  the  measured  rates  by  a  factor  of 
more  than  two.  The  application  of  Miner's  rule 

r  ^ 

resulted  in  values  of  the  order  of  .650  for 

the  constant  maximum  load  programs  and  of  .850 
for  the  variable  maximum  load  programs. 

It  is  the  utilization  of  the  facts  just  summarized,  which  were  obtained 
in  the  investigation  of  McMillan  and  Pelloux  loc.  cit.,  without  accepting  their 
conclusions,  that  we  undertake  in  the  following  pages. 
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SECTION  III 


THE  PROGRAMMED  LOAD  IN  THEORY 

We  shall  introduce  and  utilize  the  same  notation  given  in  [  5  ]•  Let 
£  be  a  subset  of  the  real-valued  piecewise  linear  functions  on  the  real  line, 
each  element  of  which  is  identically  zero  except  for  some  interval  (0,m) , 
where  m  is  a  positive  integer.  We  call  \  e  £  a  load  spectrum  whenever  the 
salient  points  of  X  occur  only  at  the  positive  integers,  the  slope  of  adja¬ 
cent  straight  line  segments  alternates  in  sign,  i.e.,  if  we  let  primes  denote 
differentiation, 

-sgn  X'(x)  =  sgn  X'(x+1)  for  all  non-integral  x>  0  , 

and  the  function  is  continuous  except  possibly  at  the  end  points  of  its  inter¬ 
val  of  support. 

The  least  positive  integer  m  such  that  x  >  m  implies  X(x)  =  0 
is  called  the  length  of  the  spectrum.  Each  unit  interval  across  which 
X'(x)  4  0  is  called  a  load  fluctuation.  The  length  is  merely  the  number  of 
load  fluctuations.  Two  successive  fluctuations,  on  one  of  which  the  load  is 
increasing,  is  an  oscillation.  Note  that  we  have  taken  the  load  oscillations 
to  be  of  the  same  duration,  i.e.,  the  frequency  of  such  oscillations  is  the 
same.  In  actual  fact  this  may  not  be  the  case  in  practice,  but  we  shall  make 
this  assumption  because  of  the  observed  invariance  of  fatigue  life  to  frequency 
in  many  cases  of  fatigue  cycling.  That  is  to  say,  we  shall  not  consider  cycles 
of  such  high  frequency  that  thermal  effects  are  encountered  or  of  such  low 
frequency  that  environmental  changes  and  age  affect  the  fatigue  life.  A  Jus¬ 
tification  of  this  assumption  has  been  made  before  by  Schijve. 
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In  order  to  make  some  of  the  following  precise  we  require  the  defini¬ 
tion  of  continuation  of  a  load  function.  Consider  Xj>Xg  €  £  with  the 

length  of  Xj  .  We  define  the  binary  operation  X^  *  Xg  on  £  by 

X1  *  Xg(x)  =  +  XgCx-mj)  whenever  X^mj)  =  Xg(0)  . 

One  checks  that  X^  *  Xg  satisfies  all  the  restrictions.  In  particular 

X  *  X  =  x2*  if  X(0)  =  X(m)  . 


We  define  a  programmed  load  as  the  repetition  of  a  spectrum  through 
its  indefinite  continuation.  It  is  thus  a  periodic  continuous  real-valued 
piecewise  linear  function  on  the  real  line  of  infinite  length.  Let  X  e  £ 
be  given  with  length  m  >  0 ;  then  X  is  the  programmed  load  with  spectrum 
X  whenever  X  =  Xj  *  Xg  *  •••  ,  where 

Xj  =  X  and  Xj+1(x)  =  Xj(x-m)  for  J  =  1,2 . 


Thus  Xj  represents  the  jth  cycle  of  the  programed  spectrum  X  .  Let 
Xjj.j  be  the  partLal  repetition  (i.e.,  the  load  history)  up  through  the  ith 
load  fluctuation  of  the  (j+1)  cycle;  that  is,  Xj+j  is  the  restriction 
of  XJ+1  in  domain  to  the  interval  (jm,  jrofi)  for  j  =  0,1,...  . 

We  now  make  the  same  probabilistic  assumptions  as  in  [  5  )>  namely, 

1)  For  given  X  e  £  the  ith  incremental  crack  extension 
during  the  last  fluctuation  of  the  loading  history  X* 
is  a  nonnegative  random  variable  Z*(X*)  depending 
only  upon  X*  .  The  random  variables  Z*(X*)  for  i  >  1 
are  statistically  independent. 
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This  assumption  implies  the  statistical  independence  of  the  crack  ex¬ 
tensions  in  each  fluctuation  not  only  from  each  other  but  from  the  total  crack 
length  as  well.  Of  course  we  do  not  preclude  functional  dependence  between 
successive  distributions,  given  the  programmed  load.  In  this  manner  the  depen¬ 
dence  upon  the  order  of  loads  is  retained. 

We  assume  as  in  [  5  ]  that 

2)  the  incremental  growth  random  variable  Z(X)  for  any 
loading  history  X  has  a  survival  distribution  R(*:X) 
which  satisfies  the  inequality,  for  x  >  0  , 


R(x:X)  I  R(t:X)dt  >  |  R(t+x:X)dt 


‘0 


‘0 


For  a  more  nearly  complete  discussion  and  intuitive  interpretation  of  this  see  [5]. 

We  now  make  the  critical  assumption,  which  replaces  the  original  inter¬ 
pretation  of  Miner's  rule  now  considered  to  be  incorrect,  namely,  that  the  load 
oscillations  may  be  permuted  in  their  order  in  any  cycle  without  altering  the 
resulting  fatigue  damage;  see  [  5  ].  This  issumption  allows  the  resolution  of 
a  load  history  into  an  equivalent  oscillation  in  terms  of  fatigue  damage,  i.e., 

3)  There  exists  a  finite  set  of  loading  oscillations,  say 

n  =  {u>0,...,uuK} 

such  that  for  any  programmed  loading  history  X*  there 
exists  an  equivalent  ujj  e  ft  for  which,  in  distribu¬ 
tion, 


Z(\l)  =  Z(u>j*)  . 

i  i* 

The  notation  will  be  X  ~  uj^ 
Consider  a  given  X  e  X  of  length  m  .  Then 
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Yj(X)  =  i  z}(X4) 

i=l 

is  the  random  crack  extension  under  the  jth  cyclic  repetition  of  X  ,  where 
different  affixes  i  and  J  on  Z*  indicate  independent  replications  of  the 
corresponding  random  variables. 

Now  we  define  the  total  crack  length  at  the  end  of  n  cyclic  repeti¬ 
tions  of  the  load  spectrun  X  as 


S„U)  =  l  Yj(X) 

J=1 

and 

EYj(X)  =  n(X)  ,  var[Yj(X)]  =  02(X)  for  j  =-  1,2 . 


It  follows  from  the  independence  that 

4n(X)  -  ESn(X)  =  nn(X)  ,  a2(X)  -  Var[Yj(X)]  =  n2a2(X) 


are  the  expected  length  and  the  variance  of  the  crack  at  the  nth  cyclic 
repetition  of  the  load  spectrum  X  . 

Let  W  be  the  critical  crack  length,  the  occurrence  of  which  defines 
failure,  and  assume 

4)  the  random  variable  W  is  statistically  independent  of 

the  crack  length  S^(X)  for  each  n  =  1,2 . 

From  these  assumptions  let  us  define  the  random  number  N(X)  of 
repetitions  of  the  spectrum  X  which  occur  before  failure  by 


[ N(X)  n]  =  t  sn-1(X)  <  W  ,  Sn(X)  >  W]  . 
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Then,  denoting  by 


EZ^Wj)  =  n(u)j)  ,  Vj  =  EN(tt)j)  for  j  =  1,*  •  ,K  , 


the  information  obtained  from  the  Wohler  diagram  of  the  distinguished  set  0 
of  loading  oscillations,  we  obtain,  by  the  arguments  in  f  5]  and  [12],  the  re¬ 
sults  which  are  a  generalization  of  Miner's  rule,  namely, 

*  n  (X) 

EN(X)  =  1/  ) 

A  J 


where 


nj(X)  =  ^  {X1  ~  (Uj*} 

i>l 

and  {«}  is  the  indicator  of  any  relation  n,  being  one  if  true  and  zero 
otherwise.  In  words  n^X)  is  the  number  of  cycles  in  the  spectrum  X 
occurring  in  their  respective  order,  which  are  equivalent  in  damage  (or  crack 
extension)  with  each  cycle  of  the  jth  load  in  continuous  repetition. 

As  a  matter  of  historical  comparison  we  recall  the  following: 

The  Miner-Pa lmgren  rule  said  that  if  K  different  load  cycles  are 
considered  and  equals  the  number  of  cycles  to  failure  under  repetition  of 

the  ith  load  cycle,  then  a  loading  spectrum  which  contains  n^  appT ications 
of  the  ith  cycle  for  i  =  1,...,h  ,  can  be  repeated  v  times  until  failure 
where 


v  = 


H 

V 

L 

i=i 


n. 
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The  onginal  interpretation  was  that  this  was  deterministic. 

A  main  purpose  of  the  subsequent  investigation  is  to  obtain  informa¬ 


tion  so  that  nj(\)  can  be  obtained  in  closed  form.  The  relationship  of 
these  concepts  to  the  distribution  obtained  in  | 12  ]  and  [ 13 )  is  readily  ap¬ 


parent  but  its  study  will  not  be  undertaken  in  this  paper. 
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SECTION  IV 


THE  THEORETICAL  BEHAVIOR  OF  FATIGUE  CRACK  GROWTH 

Let  us  consider  the  behavior  of  a  fixed  load  spectrum  X  of  length  m  , 

i.e.,  let  there  exist  m  oscillations  X^  for  1  l,«*»,m  such  that 

X  X,  *  X*-  *  •  •  •  *  X  and  the  resulting  programmed  load  X  which  has  been 
1  £  m  — 

defined  in  a  preceding  section.  After  sustaining  k  oscillations  of  the 
load  (k  1,2,*»*^  we  let  u^  be  the  expected  crack  length  and  Au^  u^  -  u^  ^ 

for  k  1  with  u^  0  .  We  set,  for  all  integers  k  , 

ak  X^k+^  *  bk  X(k)  * 


We  graph  one  such  typical  spectrum  of  m  oscillations. 


Thus  the  a^'s  represent  the  maximum  loads  and  the  b^'s  represent  the 
minimum  loads  during  the  kth  oscillation. 


✓ 
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We  now  make  the  assumption:  there  exist  two  nonnegative  functions 
and  fg  ,  each  a  f  met  ion  of  two  variables,  and  a  decreasing  sequence  of 
non- negative  numbers  ^  which  can  be  combined  to  give 

k 

4uk  [fl(Vbk*  '  L  “jfa(ak-j’  bk-j)l+  (1) 

i  1 

.L 

where  x  max(x,Ol  for  any  real  x  ,  and  without  loss  of  generality  we  take 

00 

z.  V  1  *  (2) 

j  1 

This  formulation  is  analogous  to  the  metallurgical  pronouncement  that 
fatigue  is  a  race  between  work  hardening  of  the  metal  and  the  progressive  re¬ 
laxation  of  the  yield  stress.  The  first  term  f^(a^,b^)  represents  the 
propensity  for  propagation  imposed  by  the  kth  oscillation,  which  we  assume 
can  be  expressed  as  a  function  of  the  maximum  and  minimum  stresses  of  that 
oscillation.  The  second  term  represents  the  cumulative  effects  of  work  harden¬ 
ing  or  resistance  to  crack  propagation  done  in  the  preceding  cycles.  The 
question  is  whether  such  a  simple  mathematical  model  as  presented  in  Eq(l')  can 
explain  all  of  the  various  phenomena  encountered  in  the  load  order  reactions 
of  crack  growth  which  have  been  discovered  empirically;  see  [8]  and  [  9 ] - 

The  conditions  under  which  crack  growth  was  studied  in  [8]  were  steady- 
state  conditions;  namely,  for  any  integer  k 

ac 

4uk  [W  •  L  “jf2(,'k-j))+  •  (3) 

j  i 

Let  us  now  turn  to  the  mathematical  formulation  of  the  four  empirical 
phenomena  discussed  in  [8],  to  wit,  crack  deceleration,  crack  arrest,  crack 
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acceleration  and  crack  jump.  To  dc  this  we  consider  the  behavior  of  crack 
growth  which  has  been  experimentally  determined  in  a  specified  aluminum  alloy 
(namely,  2024- T3)  under  certain  programmed  loads.  We  now  make 

Remark  1:  If  we  set  w.  0  for  i  <  0  and  define 

l  - 

ao 

l  VjnH-i  ^°r  3n^  *ntc8er  i 
j-0 

then,  if  i  <  0  and  r  is  the  least  positive  integer  such  that  rm  +  i  >  0  , 
it  follows  that  ?.  -  ?  ,  .  and 

m 

*  >  %0  >  •••  >  §  and  >  f.  1  for  any  k  l,*-*,m  .  (4) 

i-  1 

We  now  state 

Lenma  1:  The  equation  for  incremental  crack  extension  under  steady-state  con¬ 
ditions,  given  in  Eq(3),  may  be  rewritten  as 

m-1 

Auk  =  '  L  ^k-jVV^  f°r  k"  •  (5) 

J=0 

Proof:  From  Eq(3)  we  note  that 

oo  co  m 

L  wJf2^Xk-j*  =  L  't  WjnH-iVXk-jm-i)  * 

j=l  J=0  i=l 

Now  using  the  periodicity  of  X  ,  i.e.,  Xj^^  =  Xj^  for  any  integer  j  ,  we 
have  the  right  hand  side  above 

m  oo  a 

=  L  L  Wjmfif2^XnH-k-i^  "  L  ^if?^XnH-k-i^ 

i=l  J=0  i=l 

which  by  change  of  index  is  equal  to  Eq(5).  || 


15 


We  examine  first  the  programmed  loads  of  constant  amplitude  which  are 
elements  of  fi  .  Fix  co  €  fi  ;  then 

u)q  -  for  k  =  +1 ,  +2, •  •  • 


Figure  2 

Here  aQ  a^  ,  bQ  =  b^  for  k  =  +1,  +2,***  . 

From  Remark  1  and  Lemma  1  we  have  directly 


Following  our  notational  convention  f ^(uJq)  -  f^(aQ,bQ)  where  8q 
is  the  maximum  load  and  b^  is  the  minimum  load  during  load  fluctuations,  i.e., 

w0(0)  =  bQ  ,  u,0(|)  =  a0  , 
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we  assume  the  following  monotone  behavior  of  the  functions  fj  and  fg  for 
fixed  aQ  ,  bQ  >  0 

fl(x,»0)  and  fg^.bg)  are  monotone  Increasing  for  x  >  b^  , 

while 

fl(a  ,y)  and  fg(aQ,y)  are  monotone  decreasing  for  0<y<a^  . 

In  order  to  avoid  difficulties  we  shall  presume  that  the  maximum  load 
a^  is  well  under  the  yield  stress  of  the  material.  In  addition,  if  a^  =  b^  , 
there  will  be  no  load  fluctuations  (and  hence  no  fatigue),  and  as  a  consequence 
we  must  have 

1  >  ~r~7 - r  for  all  admissible  a  . 

“  fg(a,a) 

We  next  consider  the  programmed  load  which  contains  one  peak  overload. 


Figure  3 


Here  X^  =  Xj  for  i 


2,***,m-l  but  Xq  is  different,  namely, 


°o  -  ai  • 


b0^bl 


where  at  least  one  of  the  inequalities  must  be  strict;  otherwise,  Xq  =  X^  . 

It  is  known  from  [6]  that,  the  phenomena  of  crack  deceleration  and/or 
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crack  arrest  can  occur  following  such  a  peak  load.  We  now  determine  the 
behavior  of  our  model  under  a  programmed  load  of  the  type  given  in  Figure  3. 


Theorem  2:  Followiug  a  peak  overload  on  the  last  oscillation  of  a  cycle  of 

length  m,  the  crack  will  arrest  for  k<m-  1  oscillations  if  and  only  if 


f^Xj)  -  fgUj) 

W  •  w 


S  5i 


for  i  =  1 ,  •  •  •  ,k 


(8) 


and  a  necessary  and  sufficient  condition  that  the  crack  will  ad¬ 
vance  subsequently  in  the  cycle  is  that  the  reverse  inequality 
holds  for  i  k+1 , • • •  ,m-l  . 

Proof:  For  any  j  l,'**,m  we  know 

m- 1 

Auj  =  [  V  V  ‘  L  ?j-if2(Xi^ J  ' 

i-0 

Since  X^  -  Xj  for  i  l,***,m-l,  we  have 

m-1 

A“j  ffj(Xj)  -  \  Sj.jfgtXj)  -  5jf2(Xj)]+  • 
i-l 

By  adding  and  subtracting  ?j*-2^Xl^  an£*  us*n8  E<l(*0  we  obtain 

Auj  [ f1(X1)  -  f2(Xx)  -  ?j(f2(X0)  -  f2(X1))]+  .  (9) 

Again  using  Eq(U)  we  let  k  fee  the  largest  value  of  j  such  that  Eq(8) 
holds,  and  the  result  follows  by  the  monotonicity  of  .  |] 

W«*  comment  that  from  the  above  result  and  the  postulated  monotonicity 
of  f2(.,b0)  there  will  be  a  range  of  values  for  which  crack  arrest  will 
occur  and  another  range  of  values  where  the  crack  advance  will  only  be  re¬ 
duced.  We  have  the  immediate 


18 


1  2  3  4  5 

Figure  4 


If  we  set 


pj  =  Sj  -  Sj+!  for  J  =  1,  •  *  *  ,o-l  ,  Pm  =  1  +  5„  -  ?!  *  (10) 

m 

then  from  Eq(4)  we  know  o,>0  for  i=l,.  •  •  am  and  )  p  =  I  and  p  =p 

1  L  j  -  j  m- 

1 
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Under  a  progranmed  load  of  the  type  seen  in  Figure  4  we  can  state  a 
result  about  crack  acceleration  in 


Theorem  3:  If  the  maximum  stress  and  stress  variation  in  each  oscillation  of 

a  cycle  are  such  that  for  some  k  =  l,**-,m-l  we  have  Au^  >  0  and 
m 

£2^Xk- ^  £2^Xk'  (H) 

j=l 

where  for  j  =  l,**-,m  are  defined  in  Eq(  10) ,  t"»tn  the  crack 

growth  accelerates,  i.e., 

*Vl  >  A“k  •  <»« 


Proof: 

Consider 

the  condition  under 

which  Auk+1  >  Au^  ,  i.e. , 

m-1 

m-1 

fi<W 

L  ’k+l-j^V 

>  W  ■  L  ■ 

J=o 

JO 

Clearly 

W 

is  increasing  for  1 

c  l,--',m  .  Now  we  see  that 

ra-1 

m- 1 

1 

CVJ 

•n 

i 

L  ^k+l-jV  V  -  0 

JO 

J=o 

if  and  only  if 

m- 1 

L  °k-jVV  -  £2'Xk; 

.1-0 

by  adding  fgCX^)  to  sides  and  checking  the  definition  in  Eq(  10) .  Thus 

Eq(ll)  is  a  sufficient  condition  for  Eq(12).  || 

One  notes  that  Eq(ll)  would  be  satisfied  if  the  maximum  stress  and  the 
stress  variation  increased  sufficiently  in  each  oscillation. 
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SECTION  V 


A  PARAMETRIC  MODEL  FOR  CRACK  GROWTH 


Let  us  now  postulate  the  form  of  the  two  functions  which  govern  fatigue 
crack  growth;  viz.,  for  k  =  1,2  , 


^k  ^k 

f  (a,b)  =  A.  a  (a  -  b)  for  a  >  b  >  0 


(13) 


where  A^  ,  ,  3^  for  1  =  1 ,2  are  unknown  parameters,  assumed  to  be 

positive,  which  must  yet  be  determined.  It  is  clear  that  these  functions 
both  satisfy  the  monotonicity  requirements  specified  in  the  preceding 
section. 

From  Theorem  1  we  know  that  for  a  programmed  load  of  constant  maximum 
stress  a  and  minimum  stress  b  the  satisfaction  of  the  condition 

3,-0,. 

a  A  ^(a  -  b) 

is  a  necessary  and  sufficient  condition  that  the  crack  advance.  We  thus 
assume  that 


1  < 


f1(a,b)  Al  Oj-Og  .  N K1’K2 


f^TaTbl  A, 


0  <  Aj  <  Ag  ,  >  Q!g  >  0  ,  3j  >  Pg  >  0  . 


W 


We  note  that  as  a  consequence  of  the  formulation  of  Eq(l3)  the  "endurance 
limit",  i.e.,  a  value  of  the  maximum  stress  below  which  fatigue  will  never 
occur  (we  say  "never"  because  in  this  model  we  are  not  considering  any 
stochastic  variation  in  the  damage  caused  by  each  oscillation  of  the  load), 
will  be  altered  if  one  takes  different  values  of  the  minimum  stress  unless 
we  assume  3,  -  3-  .  One  of  the  consequences  of  the  basic  stochastic  nature 

i  C 

of  the  fatigue  process  is  that  such  questions  us  whether  or  not  some  parameter 
is  zero  or  just  very  small  cannot  be  resolved  by  physical  measurement. 
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Subsequently  we  propose  to  '  the  simplest  mathematical  assumptions 


possible  which  are  consistent  with  the  facts. 

ihe  crack  growth  rate  per  cycle,  under  constant  amplitude  loading, 
assuming  it  takes  place,  is 

“*1  ^2  ^2 
Au  =  f^a.b)  -  fg(a,b)  =  A^a  (a-b)  -  Aga  (a-b) 


The  claim  of  reference  [  11  ]  is  that  for  a  fixed  maximum  stress  the 
logarithm  of  the  crack  growth  rate  per  cycle  is  nearly  a  linear  function 
of  the  logarithm  of  the  stress  variation.  However,  the  evidence  also  shows 
that  this  slope  increases  as  the  variation  in  stress  decreases  (i.e.,  the 
slope  is  a  decreasing  function  of  the  variation  in  stress);  see  p.  33 
loc.  cit. 

Letting  x  =  <Ln(a  -  b)  with  a  fixed,  we  find  as  a  function  of  x 


3.x 

Au  =  Aje  [  1  -  Ke  01  ] 


(15) 


for  some  constants 


a. 


Aj  =  Aa 


A2 


K  =  A"  a 
A1 


e  -  Pj  -  P2  >  0  . 


(16) 


We  set 


d  Au 


•P(x) 


and  we  check  that  p'(x)  <  0  . 

We  have  verified  that  the  model  Eq(13)  behaves  in  this  regard  as  does 
the  actual  fatigue  crack  itself  as  recorded  in  [11].  Now  we  verify  another 
experimentally  determined  condition  concerning  how  the  growth  rate  behaves 
as  a  function  of  the  maximum  stress  when  the  stress  variation  is  held  constant. 
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It  is  stat'd  (p.  33  loc.  cit.)  that  the  logarithm  of  the  crack  growth  rate 


is  again  at  'mately  a  linear  function  of  the  maximum  stiess  but  that  the 
slope  is  a  decreasing  unction  of  this  maximum  stress. 

Letting  y  =  lr\  a  with  fixed  difference  a  -  b  ,  we  have 


Au  =  A'j'e 


V 


[1 


K’e'6y] 


(17) 


where 

f^l  Ag  ^2"^1 

A'l  =  Ajl3  “  b)  1  >  0  ,  K'  =  ^(a-b)*  >0, 

and  6  =  0:^  ■  a2  —  ®  are  constants.  By  comparison  with  the  preceding  results 
we  see  that  the  requisite  behavior  of  the  growth  rate  follows  also  in  this  case. 

Where  we  disagree  with  the  conclusion  of  McMillan  and  Pelloux,  p.  57 
#3  in  [11]  is  in  taking  the  crack  growth  Au  to  be  of  the  same  parameteric 
form  as  fj  and  f^  .  Not  only  does  this  form  seem  unsubstantiated,  but  it 
does  not  allow  for  crack  retardation. 

The  similarity  of  these  two  cases,  Eq(l5)  and  Eq(l6),  prompts  us  to 
propose  the  same  method  for  the  determination  of  the  unknown  parameter  for  the 
case  of  constant  amplitude  loading. 


Theorem  2:  A  necessary  and  sufficient  condition  that  there  exist  parameters 
e  >  0  ,  K  >  0  such  that  the  convex  function  q  defined  by 
Ke 


q(x)  = 


ex 

e  -  K 


for  c  <  x  <  d 


(18) 


satisfies  the  conditions  q(c)  =  C  >  q(d)  =  D  is  that 

.  1  1 
d  -  C  =  D  -  C  • 


(19) 


In  that  case  e  is  the  solution  of  the  equation  g(y)  =  0  ,  where 
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g(y)  -  [1  +i]De'’(d-c>  -  D  -  y  , 


(20) 


and  in  terms  of  e 


K  = 


Ce 


ec 


C  +  e 


(21) 


Proof:  The  two  conditions  q(c)  =  C  ,  q(d)  =  D  yield  the  two  equations 
C(e€C  -  K)  =  Ke  ,  -  K)  =  Ke  . 


Solving  for  K  in  the  first  equation  yields  Eq(21).  Using  this  to  eliminate 
K  in  the  second  yields  g{  e)  =  0  with  g  defined  in  Eq(20) .  Now  g  is  a 
convex  function  such  that  g(0)  =  0  ,  g(°°)  =  00  .  For  e  >  0  we  must  have 
g'(0)  <  0  ,  which  is  the  condition  of  Eq(l9).  || 


Lemma  2: 


If  is  given,  more  accurate  approximations  can  be  found  by 

Newton- Raphson  iteration  for  n  =  0,1, •••  : 


-e  (d-c) 

e  (d  -  c)[  1  +  e  / C]  -  1  +  e 
n  n 

(d - c) [ 1  +  en/C]  +  1/c  -  ^e  Cn^d 


(22) 


and  c  -*  e  as  n  -*  00  . 

”  *(«n) 

Proof:  This  follows  directly  from  e  =  e  -  -77 — N  and  the  nature  of  g  .  II 

7  n+l  n  8  (en) 

From  the  extant  data  we  now  seek  to  establish  values  of  the  parameters 
which  determine  both  f^  and  fg  .  We  first  attempt  to  analyse  the  results  of 
programmed  spectra  with  constant  maximum  loads  and  variable  amplitude.  Unfortu¬ 
nately  measurements  of  incremental  crack  advance  for  programed  spectra  of  only 
constant  amplitude  loading  are  not  available  in  [11]  and  comparable  experimental 
results  of  sufficient  accuracy  are  not  readily  available  elsewhere.  Moreover, 
in  [11]  the  crack  growth  rates  for  constant  maximum  loads  with  variable  amplitude 
but  with  load  order  interchanged  were  not  separated,  and  the  data  on  the 


24 


slope  of  the  logarithm  of  the  crack  advance  rate  were  put  together.  This 
Ct.  founding  of  the  data,  we  believe,  does  not  interfere  in  our  determination 
of  the  parameters.  This  is  one  of  the  more  interesting,  or  fortuitous,  re¬ 
sults  of  this  investigation,  especially  since  the  inference  from  the  model 
agrees  with  the  empirical  conclusions  of  other  investigators. 
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SECTION  VI 


THE  DETERMINATION  OF  PARAMETERS  USING  SPECTRA  OF 
CONSTANT  MAXIMUM  STRESS  BUT  VARIABLE  AMPLITUDE 

We  first  quote  a  statement  of  McMillan  and  Pelloux  [11]  p.  33.  1.8 
et  seq.  concerning  spectra  of  the  type  under  discussion.  "There  was  a  certain 
scatter  in  the  data  but  it  showed  that  for  constant  S^^  (maximum  stress) 
the  following  relative  crack  rate  equation  could  be  written 

rate  =  constant  (AS)n  (23) 

(where  AS  is  the  stress  variation)  with  n  =  1  for  AS  between  eight  and 
twelve  thousand  psi  and  increasing  from  two  to  four  for  AS  lower  than  eight 
thousand  psi." 

By  comparing  Eq(l5)  with  Eq(23)  we  see,  in  our  terminology,  we  wish 
to  choose  the  free  parameters  so  that,  letting  x  denote  the  <tn-stress 
range,  we  have 

p(x)  >  1  for  <ln(8xl03)  <  x  <  ^(12x10^) 

but  moreover  so  that  p(x)  increases  to  between  2  and  4  whenever 
o 

x  <  £n(8  x  10J)  .  In  order  to  utilize  Lemma  2  we  set 

d  =  <tn(12xl03)  c  =  £n(6xl03) 

P(d)  -  3/2  P(c)  =  3  (24) 

D  3/2  -  0=5-^. 

Unfortunately  the  scatter  in  the  data  will  prohibit  an  easy  and  exact 
determination  of  the  disposable  parameters.  (This  is  indicated  in  the  para¬ 
graph  following  the  one  quoted.)  Nevertheless,  so  that  we  retain  the  mono¬ 
tonicity  which  p( • )  should  exhibit  over  its  range* our  latitude  in  the 
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selection  of  P^  must  be  restricted  in  order  for  its  behavior  to  remain  con¬ 
sistent  with  the  facts. 

We  first  must  check  that  Eq(19)  is  satisfied,  namely, 

(3-PjMf-Pj)  <  2^2  =  2-l64(*  • 


Analysis  reveals  that  this  inequality  is  satisfied  for  0.6l8  <  Pj  <  3*88  . 
In  order  that  D  >  0  we  must  keep  0.6l8  <  Pj  <  1.5  •  It  follow^  since 
c  =  Pi  *  and  0g  >  0  ,  that  e  <  Pj  .  Since  g  is  convex,  a  necessary 
and  sufficient  condition  for  this  is  that  g(P^)  >  0  .  Substituting  Eq(24) 
into  Eq(20),  we  see  that 


g(P^)  >0  if  and  only  if 


»jtn2  3  -  Pj 
e  -  F237  ‘ 


One  sees  that  this  is  true  for  0  <  Pj  <  1  .  Thus  the  permissible  range  of 
P^  is  only 


0.62  <  Pj  <  1  . 


(25) 


We  pick 

Pj  =  1  D  =  \  ,  C  =  2  . 

Utilizing  Lemma  2  to  calculate  €  ,  we  first  guess  =  1  and, 
to  our  surprise,  calculate 

=  3«n  i)-it  =  , 

1  3(<tn  2)  +  1  -  4e-*n  2 

Therefore  we  have 


e  =  1  ,  p2  =  0  .  (26) 

From  the  determination  just  made  we  have 
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Lemma  3:  For  maximum  stress  a  and  minimum  stress  b  the  functions 
and  fg  extending  and  resisting  crack  growth  are 

“l  a2 

f1(a,b)  =  Aja  (a-b)  ,  fg(a,b)  =  k^a 

where  0:^  ,  are  yet  unknown  parameters. 

Thus  we  have  the  following  conclusion  on  the  behavior  of  crack  growth 
which  we  state  in 


Theorem  3:  Under  a  spectrum  X  of  constant  maximum  stress  and  variable  stress 
amplitude,  i.e., 

X(k+|)  =  a 

k  -  1 , •  •  •  ,m 

X(k)  =  bk 

the  crrck  growth  for  the  kth  oscillation  is  independent  of 
,?m  and  is  given  by 

°2  + 

Auk  =  ^Ala  ^a'bk^  ’  A2a  ^  k  m  . 


Proof:  From  the  general  result  on  steady  state  crack  advance  in  Lenina  1  we  have 


ra- 1 

v 


IW  ■  f2<V  l  W 

3=o 

but  by  Lemma  3  we  know  t0  depends  only  upon  the  maximum  stress  and  not  upon 

*  m-1 

the  stress  range.  But  for  all  k  =  l,***,m  we  know  E  §.  =  1  and  we 

J=0  k'J 

have  the  result.  I 


A  conclusion  of  this  theorem  is  that  any  spectrum  with  a  constant 
maximum  stress  at  each  oscillation  but  variable  minimum  stress  would  yield 
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the  same  total  crack  advance  per  spectrum  independently  of  the  order  in  which 
the  loads  are  applied. 

This  mathematical  result  coincides  with  the  empirical  conclusion  of 
p.  30,  1.  14  et  seq.  [11]: 

"This  shows  that  in  the  tests  at  constant  maximum  load  the  sequence 
of  load  application  did  not  measurably  influence  the  overall  crack  growth 
rates.  This  was  confirmed  by  the  fracture  analysis  with  the  electron  micro¬ 
scope: 

1.  There  was  no  marked  crack  front  advance  at  the  change  of  load 
amplitude  between  different  spectra. 

2.  After  a  change  of  load  amplitude,  the  spacing  of  the  striations 
for  the  next  load  amplitude  sequence  reached  a  stable  and  uniform  value  on 
the  first  cycle  of  the  new  load  amplitude  sequence." 

We  conclude  this  section  with  the  following  calculation:  Since 
c  =  ln(6  x  103)  for  this  series  of  experimental  loading  programs,  PI,  P2, 

P3,  P4,  we  substitute  the  known  parameters  into  Eq(21)  to  obtain 

„  Ceec  2  tn(6xlo3)  ,,  lr3 
K  =  =  3  e  =  4x  1(T  .  (27) 

From  the  definition  of  Eq(l6)  and  the  fact  that  the  constant  value  of 
the  maximum  load  was  a  =  12  x  1C  ,  we  have  a  determination  of  the  ratio 

^  =  k(12x  103)6  .  (28) 

A1 

Once  a  value  of  6  is  determined  this  ratio  cau  be  calculated. 
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SECTION  VII 


THE  DETERMINATION  OF  PARAMETERS  USING  SPECTRA  OF 
CONSTANT  STRESS  VARIATION  BUT  VARIABLE  MAXIMUM  STRESS 

We  quote  a  conclusion  of  McMillan  and  Pelloux  [  11]  p.  39>  1*  8 
et  seq.  with  regard  to  programs  P6,  P7 ,  P8,  P9- 

"The  average  microscopic  rates  of  crack  growth  per  load  cycle  for 
Spectra  A  and  B  of  Programs  P6,  P?,  P8,  and  P9  were  measured  within 
each  program  and  compared  with  the  larger  growth  rate  (taken  as  unity)  cor¬ 
responding  to  S  max  14,000  psi.  The  effect  of  the  maximum  load  amplitude 
on  the  relative  crack  growth  rate  per  cycle  is  shown  in  Figure  6  .  The 
scatter  was  quite  large,  but  a  relationship  of  the  type: 

rate  =  constant  (S  max)m 

could  be  written  to  represent  relative  crack  growth  rate.  At  a  given  crack 
length  and  for  a  constant  maximum  load,  m  varied  from  two  to  three  at  large 
S  max  to  three  to  five  at  smaller  S  max.  In  the  same  figure  the  relative 
microscopic  rates  confirm  the  fact  that  P9  had  a  larger  growth  rate  than 
P8.  This  can  be  explained  by  assuming  a  crack  growth  retardation  in  Spectrum 
B  of  P8  because  Spectrum  B  followed  the  higher  load  level  of  Spectrum  A  . 
The  crack  growth  retardation  could  be  due  to  a  work  softening  of  the  plastic 
zone  by  Spectrum  B  and  a  consequent  lower  crack  growth  rate." 

We  now  attempt  to  utilize  the  same  procedure  employed  in  the  preceding 
case  to  determine  the  unknown  parameters.  Unfortunately,  there  are  no  data 
on  crack  growth  rates  for  spectra  each  of  the  same  constant  stress  variation 
but  with  different  and  constant  maximum  stress.  The  fractographic  difficulty 
with  analyzing  such  a  spectrum  would  be  that  without  some  measurable  stress 
variation  in  each  program  the  appropriate  section  of  the  fracture  face  could 
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not  be  made  to  correspond  with  the  imposed  load.  However,  since  only  sections 
of  spectra  of  constant  stress  variation  could  be  identified  with  the  crack 
growth  measured,  there  is  certainly  an  influence  due  to  the  preceding  blocks 
of  higher  or  equal  maximum  stress. 

Consider 

m-1  ^ 

%  =  lW  '  L  VjWJ  • 

J-o 

Since,  by  definition,  a^  -  b^  is  constant  for  k  =  l,-**,m  ,  we  set 

A'i  =  Ai(ak' V  *  K'  =  VA'i  ; 


then  utilising  the  results  of  Theorem  3  we  have 

m-1 


a 


Au,  = 
k 


ir 


Al\  L1  -  L  5k-j"  -j 
j=0 


or  -a  + 

K'a , 


Let  a  denote  the  constant  value  of  the  maximum  stress  within  the 
block  in  the  spectrum  which  is  of  interest.  By  examining  the  programmed 
loads  P6,  P7 ,  P8,  P9  we  see 


m- 1 

\  -  -6  -6 
L  5k-jaj  i  a 
J=0 


where 


6  - 


a. 


It  follows  that 

ai  _* 

Au^  <  A'jS  [  1  -  K'a  ]  . 

Thus  the  true  crack  growth  rate  within  the  block  would  be  less  than  that  in 
a  corresponding  constant  amplitude  test  with  the  same  maximum  stress.  To 
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proceed  in  the  same  manner  as  before  we  write  . 


Au  =  A'^a  *[1  -  K'a'6] 


and  letting  y  =  in  a  we  have 


9  Au  _  .  K '  6  _  ^ ,  \ 

■  “i +  7777: ' p  (y)  • 


In  order  to  utilize  Lemma  2,  we  attempt  to  compensate  for  the  reduction 
in  the  rate  as  demonstrated  above.  Accordingly  we  interpret  the  data  to  adjust 


by  setting 


c  =  -tn(lC>  ) 


d  4n(  1.4x10) 


P*(c)  =  5 


p*(d)  =  3 


c  =  5  -  a. 


d  =  3  -  a. 


We  now  check  that  Eq( 19)  is  satisfied,  namely, 

4n(!.4)  <  -  5  .  ^  . 

Analysis  shows  this  inequality  satisfied  for  1.355  <  <  6.644  but  since  we 

must  also  have  D  >  0»  we  require  Q!^  <  3  .  As  before  it  follows,  since 
6  =  ai  '  a2  and  Ctg  >  0  *  that  6  <  Q:^  .  A  necessary  and  sufficient  condition 
for  this  is  that  g(c^)  >  0  with  the  proper  parametric  values  substituted 
into  g  .  Using  Eq(20)  we  see 

a  Lnl.4  .5  -  a  . 

g(0'1)  >0  if  and  only  if  e  >  J  . 

Numerical  analysis  shows  this  to  be  true  for  0  <  ctj  <  2.2-  .  Thus 
the  combined  restriction  on  Q:^  is 


1-35+  <  otj  <  2.2- 
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We  pick 


Oj  =  2  C  =  3  »  D  =  1 


(30) 


We  now  check,  using  Lemma  2,  that,  if  we  select  a  value  6q  =  1.5  » 
we  find  6^  =  1.4966  .  Thus  we  take 


6  =  3/2  £^  =  ^-6=1/2  . 

We  recall  that  c  =  im(lO^)  and  hence 


6c 


(31) 


K*  _  Ce^_  2  6 

c+6  ~  3xl°  * 


Since  the  stress  variation  was  a^  -  =  7  x  103  psi*  we  have 


K'  =  -£(7x  103)’1 
A1 


and  thus 


^  =  K'(7x  103)  -  4.666  x  109  . 


However,  we  have  precisely  calculated,  in  Eq(27),  that 

~  =  K(  12 x  103)6  =  5-259  x  109 
A1 

the  last  equality  obtained  by  using  6  =  ^  and  the  value  of  K  from  Eq(26) . 
There  is  a  discrepancy  of  only  about  13  1>  between  these  two  determinations. 

(The  author  does  not  pretend  that  these  virtually  coincidental  values  were 
obtained  on  the  first  iteration.) 

q 

The  average  value  of  the  two  determinations  of  Ag/A^  is  4.962x10^  , 
and  thus  for  our  future  calculations  we  take  the  ratio 


2  x  10 


-10 


(32) 


33 


The  disparity  between  these  two  determinations  could  be  reduced  by  ad  hoc 
adjustment  of  the  initial  values  and/or  at  the  expense  of  the  rational  ex¬ 
ponents  within  the  model;  however  the  appeal  of  such  values  comes  not  only 
from  their  simplicity  but  also  from  the  possibility  that  the  basic  mechanism 
of  fatigue  can  be  explained  in  simple  mechanical  terms,  whatever  the  material. 
From  our  determination  of  the  exponents  we  obtain 

Lemma  4:  For  a  maximum  stress  &  and  a  minimum  stress  b  within  a  load  oscil¬ 
lation  we  have,  for  the  functions  extending  and  resisting  crack  growth, 

f^a.b)  =  a2(a-b)  ,  fg(a,b)  =  /a  . 

We  now  state 

Theorem  4:  Under  a  spectrum  X  of  length  m  with 

X(k  +  |)  =  ak  ,  X(k)  =  bk  for  k  =  l,---,m 
the  crack  growth  resulting  from  the  kth  oscillation  is  given  by 

a  m_1  + 

=  -  I  ?k.j*kJ  •  ‘33> 

J=0 

Note  that  the  appropriate  values  of  ^  have  yet  to  be  determined. 

Turning  our  attention  to  the  actual  dimensions  of  the  variables  involved 
in  the  calculations,  we  find  the  measurements  of  Au  ,  the  crack  growth  rate,  are 
expressed  in  microns  which  range  in  value  from  0.1  to  1  whenever  a  ,  the  maximum 
stress,  is  measured  in  ksi  (ksi  equals  a  thousand  psi)  varying  from  5  to  15  ,  and 
(a  -  b)  ,  the  stress  variation,  is  also  expressed  in  ksi,  varying  from  2  to  12  . 
This  mixture  of  metric  and  English  systems  results  from  the  calibration  of  the 
test  instruments  used.  However,  consistent  dimensioning  can  be  assured  by  proper 
choice  of  units  to  express  the  constants  and  Ag  . 
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We  have  thus  so  far  determined  that  under  the  conditions  extant, 
namely,  202k-T3  aluminttm  alloy  cycled  in  laboratory  air  at  room  temperature 
with  stress  being  expressed  in  tens  of  ksi  and  a  crack  length  between  0.6  and 
1.7  inches,  the  incremental  crack  advance  is  given  by 

m-1  + 

‘"k  =  M2**' V  V  •  I  WajJ  (34) 

J=0 

for  steady  state  advance  at  the  kth  oscillation. 

Note  that  the  constant  of  proportionality  and  the  "memory"  con¬ 

stants  £  (i  =  l,...,m)  must  be  determined  by  other  measurements  using 
different  data. 

A  comment  is  now  in  order.  The  mechanism  of  resistance  to  crack  growth 
results  from  the  plastic  yield  zone  which  precedes  the  crack  tip.  The  size 
and  shape  of  this  zone  are  determined  by  the  loads  of  the  spectrum  which  have 
been  previously  applied,  as  well  as  the  geometry  of  the  specimen.  It  is 
admittedly  a  simplification  to  suppose  that  the  crack  growth  resistance  can 
be  expressed  as  a  one  parameter  functi  n  of  the  number  and  magnitude  of  the 
loads  preceding.  We  feel  that  this  simple  supposition  yields  a  more  realistic 
description  of  the  phenomenon  than  any  model  proposed  heretofore. 

We  let  the  weighting  factors  be,  for  some  p  <  p  <  1  , 

Wj  =  (1  -  p)pl  1  i  =  1,2, •••  . 


Hence  for  a  spectrum  of  length  m  we  have 


r  =  '  p)?1-1 

'i  .  m 

1  -  p 


m 

and  clearly  E  ?,  =  1  . 

1  1 


i  =  1,- 


•,m 


(35) 
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SECTION  VIII 


CONCLUSIONS 


In  this  study  we  have  concentrated  on  the  problem  of  determining  a 
formula  for  the  effect  of  load  order  interactions  on  the  incremental  fatigue 
crack  growth  resulting  from  the  load  oscillations  in  a  loading  spectrum.  We 
have  utilized  only  data  on  aluminum  2024-T3  alloy. 

We  have  argued  that  the  incremental  growth,  if  of  the  form 


r  ai  '  ao 

4uk  =  LVk  K  ■  V  *2  l  5k-jai  (ivbj)  J  ■ 

j-o 

where  is  the  maximum  stress  and  is  th^  minimum  stress  for 

the  kth  load  oscillation  with  as  memory  constants,  will 

satisfy  all  the  observed  phenomena  of  crack  growth, namely,  crack  jump,  crack 
arrest,  and  crack  acceleration  and  deceleration. 

Since  there  is  no  load  order  interaction  in  the  crack  growth  data 
obtained  when  the  maximum  stress  is  held  constant,  the  determination  of  the 
parameter  values  0^  =  1  and  0g  =  0  by  the  method  chosen  would  seem  to  have 
more  assurance  than  the  determination  made  subsequently  of  and  from 

the  data  obtained  when  the  stress  variation  was  held  constant. 

This  is  because  we  must  guess  at  the  reduction  caused  through  load 
order  interactions,  i.e.,  the  net  effect  of  the  reduction  factor 

m-l  g 

E  §,  a.  with  6  =  a.  -  cr_  ,  in  choosing  the  values  of  the  slope  of  the 
j=0  K'J  J  1  * 


growth  to  be  somewhat  higher  than  the  central  values  indicated  by  the 
data.  This  has  resulted  in  the  values  =  2  ,  =  1/2  .  The  necessity 

of  this  extrapolation  is  a  consequence  of  the  type  of  spectrum  which  was 
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studied.  It  is  clear  that  those  experiments  in[ll]were  not  run  with  the 
idea  of  estimating  the  parameters  of  such  a  model  as  we  have  proposed  but 
merely  of  strengthening  the  hypothesis  that  crack  structures  could  be  directly 
related  to  the  load  oscillations  within  each  programmed  load. 

It  is  conceivable  that  subsequent  data,  more  suitable  for  estimation, 
will  reveal  that  these  chosen  values  were  in  error.  However,  since  there  are 
additional  parameters  which  are  adjusted  to  fit  the  actual  growth  data,  such 
a  discrepancy  could  only  be  shown  by  prediction  from  our  model  to  other 
situations. 

We  have  determined  the  following  formula  for  incremental  crack  growth 
for  the  kth  load  oscillation,  based  on  the  somewhat  ambiguous  data  (from 
our  point  of  view)  of  [11]: 


2  m”1 

Auk  =  A^a^b^  -  £  5k_/aj] 

JO 


k  =  1, 


where  a  and  b  are  both  expressed  in  tens  of  ksi. 

We  note  that,  if  we  had  =  3  (instead  of  2),  then  for  sinusoidal 

loading  the  incremental  crack  growth  would  be  approximate!'  jf  the  form 
4 

Au^  -  constant  a^  ,  which  would  be  consistent  with  the  so  called  fourth  power 
rule  given  by  Paris  [14], 

On  the  basis  of  the  above  formula,  if  is  the  set  of  constan 

amplitude  spectra  of  different  maximum  stress  for  which  the  mean  number  of 
cycles  until  failure  is  known,  then  for  a  given  programmed  spectrum  X  with 
history  until  the  kth  oscillation,  namely. 
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we  have  the  equivalence  relation 


.  k  k 
X  ~  (l)j 


if  and  only  if 


Auk(>.) 


Thus  the  formula 


m 

nj(X)  =  £  [Auk(X)  =  Auk(tuj)} 
lcl 

can  be  used  to  determine  the  expected  life  under  repetition  of  the  spectrum  X  . 

To  show  the  overall  agreement  of  the  fitted  crack  growth  model  given 
in  Eq(3^)  with  the  observed  behavior,  we  exhibit  in  each  of  Figures  7,  8,  9> 

10  and  11  two  calculations  (made  with  two  different  values  of  the  memory 
constant  p  as  defined  in  Eq(35))  of  the  incremental  crack  growth  resulting 
from  the  programmed  loads  displayed.  We  also  show  the  actual  fracture  surface 
of  the  aluminun  alloy  upon  which  the  load  was  imposed.  The  photographs  have 
been  reproduced  from  reference  [11]  with  the  authors'  permission. 


38 


Relative  Crack  Growth  Rates 


Figure  5 

Scatter  of  Relative  Crack  Growth  Rates  at  Constant  Maximum 
Loads  for  Programs  Pi  ,  P2  ,  P3  ,  P4 


Relative  Crack 
Growth  Rates 


Figure  6 

Scatter  of  Relative  Crack  Growth  Rates  at  Constant 
Load  Amplitude  but  Variable  Maximum  Load 


Figure  7 

Typical  fracture  surface  due  to  Program  P  2,  showing  the  uniformity 
of  striation  spacing  during  each  of  the  first  two  groups  of  oscilla¬ 
tions.  Note  the  mathematical  model  predicts  crack  arrest  during  the 
third  group  of  oscillations. 
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Figure  8 


Typical  fracture  surface  due  to  Program  P  6.  Note  the  large  crack 
jump  during  the  first  oscillation,  which  the  mathematical  model 
predicts . 
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Figure  9 

Fracture  surface  typography  due  to  Program  P  10.  Note  the  sharpness 
of  striation  ridges  indicating  a  sharp  crack  tip.  The  difficulty  of 
identifying  striations  with  actual  loads  is  apparent  of  random  loads 


Figure  10 

Typical  fracture  surface  resulting  from  Program  P  11.  Note  in 
the  photograph  the  large  striation  spacing  due  to  the  load 
amplitude  in  the  13th  oscillation,  followed  by  three  large 
striations  corresponding  to  the  next  three  oscillations.  We 
see  the  mathematical  model  agrees  with  such  behavior. 
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Figure  11 

Typical  fracture  surface  resulting  from  Program  P  12.  Note  the  large 
striation  spacing  due  to  the  load  oscillation  in  14th  cycle,  preceeded 
by  three  smaller  striations.  The  mathematical  model  agrees  with  this 
behavior. 
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LIST  OF  SYMBOLS 


Z 


a  set  of  functions  of  time,  the  values  of  which  give  the 
deflection  of  a  metallic  specimen  being  fatigued 


X  =  X.*  •  •  **x 

1  m 

X 

Z(X) 

Yj(X) 

sn(X) 

EYj  (  X)  =  n(X) 

Var[  Yj(X)  ]  --  o2(X) 
N(  X) 


the  jth  load  oscillation 

a  load  spectrum  consisting  of  m  load  oscillations 

a  programmed  load  (or  spectrum)  consisting  of  an  indefinite 
number  of  repetitions  of  X 

the  (random)  incremental  crack  growth  following  the  last 
fluctuation  in  X 

the  finite  set  of  reference  loading  oscillations 

the  total  crack  extension  resulting  from  the  jth  repetition 
of  spectrum  X 

the  total  crack  extension  at  the  end  of  n  cyclic  repeti¬ 
tions  of  spectrum  X 

the  expected  (mean)  crack  extension  per  sper trum 

the  variance  in  crack  extension  per  spectrum 

the  (random)  number  of  repetitions  of  spectrum  X  which 
can  be  repeated  until  fatigue  failure  occurs 


a^  ,  b^  the  maximum  and  minimum  values,  respectively,  of  the  stress 

during  the  kth  load  oscillation  in  the  spectrum 

f.(X.)  the  forcing  function  for  crack  extension  during  the  kth 

fluctuation 


fg(X^)  the  resisting  function  for  crack  extension  during  the  kth 

fluctuation 


w . 


the  weighting  factor  which  represents  the  "work  hardening" 
at  the  crack  tip  resulting  from  the  jth  preceding  load 


Au^X) 


the  mean  incieraental  crack  extension  during  the  kth 
fluctuation  in  the  spectrum  X 


the  cumulative  weighting  factor  for  the  ith  load  in  the 
spectrum  during  steady  state  growth 

f^(a,b)  for  i  1 ,2  the  resistance  or  forcing  function  expressed  in  terms  of  the 
maximum  stress  a  and  minimum  stress  b 
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y=(na 
x  =  -tn(a  -  b) 

p(x) 

a±  for  i =  1,2 
(3^  for  i =  1 ,2 


the  natural  logarithm  of  the  maximum  stress 

the  natural  logarithm  of  the  stress  variation 

the  slope  of  the  In- incremental  crack  growth  as  a  function 
of  x 

the  exponent  of  maximum  stress  in  f^ 
the  exponent  of  stress  variation  in  f^ 


e  = 


A 


with  or  without  affixes  represents  a  constant  associated 
with  incremental  crack  growth  behavior 


q(x)  =  p(x)  -  3X 

g(y)  =0  an  equation  the  solution  of  which  gives  the  value  of  e 

satisfying  specified  boundary  conditions 


c  ,  d  ,  C  ,  D 


constants  which  represent  the  boundary  conditions 
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